Vortex Dynamics: Quantum versus Classical Regimes 
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For many years the classical Hall-Vinen-Iordanski (HVI) equation has been used to analyse vortex 
dynamics in superfluids. Here we discuss the extension of the theory of vortex dynamics to the 
quantum regime, in which the characteristic vortex frequency is higher than the temperature. At 
the same time we justify, in the low-frequency classical regime, the use of the HVI equation, provided 
an inertial mass term and a noise fluctuation term are added to it. The crossover to the quantum 
regime is discussed, and an intuitive picture is given of the vortex dynamics, which in general is 
described by 2 equations (one for the vortex coordinate, and one for its quantum fluctuations); we 
also discuss the simple equation of motion found in the extreme quantum regime. 

PACS numbers: 



INTRODUCTION 

Very soon after the discovery of quantum vortices in 
superfiuid 4 He by Vinen [T] , an equation of motion for a 
vortex was proposed by Hall and Vinen[5]; a few years 
later Iordanskii 3J added an extra term, to produce what 
is commonly known as the Hall-Vinen-Iordanskii (HVI) 
equation of motion. The HVI equation, occasionally sup- 
plemented by an inertial mass term and by a noise fluc- 
tuation term, has been used in the last 60 years to ana- 
lyze thousands of experiments in superfluids and super- 
conductors. However it has been controversial, and for 
the last 20 years a strenuous debate has been going on 
over its validity. Key questions concern the value of the 
vortex effective mass M v (estimates range from zero to 
infinity jl]) and the vortex-quasiparticle coupling coeffi- 
cients D (T), D' (T). Indeed, Thouless et. al. [5] find 
D' Q (T) — for all T; and scattering analyses [3l |6H2] give 
various different results for D' a (T). 

In what follows we wish to shift the focus of this dis- 
cussion. In our view, the key question is how to give 
a proper quantum-mechanical description of the vortex 
dynamics. This requires two equations of motion, one 
for a vortex 'centre of mass' coordinate, and the other 
for the quantum fluctuations around this coordinate [10]. 
A key parameter in the theory is the dimensionless ra- 
tio = tel/kT, where T is the temperature, and Q the 
characteristic frequency of the vortex dynamics. In the 
classical limit, where O — >■ 0, one actually recovers the 
HVI equation with added inertial and noise terms [TO]. 
While key questions remain (notably the role of bound- 
aries and normal fluid velocity |llj . and the calculation 
of the effective mass [J]), we believe that the main points 
at issue have now been settled in the classical regime. 

However the quantum regime remains relatively un- 
explored. It possesses a number of fascinating features, 
both theoretical and experimental. The theory gives an 
illuminating picture of vortex-quasiparticle interactions, 
classical analyses of which interactions have been subtle 



and controversial, because of the long-range nature of the 
interaction, and the difficulty of accounting for vortex 're- 
coil' in the scattering. A quantum analysis immediately 
makes clear that a key feature of the interaction (which 
is properly described as an interaction between quantum 
soliton and quasiparticle field excitations in 2 + 1 dimen- 
sions) , is the distortion of the quasiparticle part of the su- 
perfiuid wave-function by the vortex 'zero mode' part of 
this wave-function. On the experimental side the results 
suggest new kinds of experiment, and show that most 
previous experiments have been in the classical regime. 

In this paper we present a more intuitive picture of the 
quantum regime for vortices in a Bose superfiuid, high- 
lighting the new features. On the theoretical side we first 
discuss the main features of a quantum description of the 
vortex, and of vortex-quasiparticle interactions. We then 
briefly summarize the resulting equations of motion for 
the vortex, and give an intuitive picture of them in real 
time (as opposed to frequency space). Almost all details 
of the calculations are eschewed - they are quite lengthy 
and will appear elsewhere [H]. We then very briefly dis- 
cuss how one might experimentally probe the quantum 
regime in a Bose superfiuid. We emphasize that this pa- 
per is about vortices in neutral Bose superfluids - the 
problem of vortex dynamics in Fermi superfluids (both 
charged and neutral) is rather different (and is, we be- 
lieve, still open). 



QUANTUM-MECHANICAL DESCRIPTION OF A 
SUPERFLUID VORTEX 

In a classical description, a superfiuid vortex is de- 
scribed by its position R„(t) (and time derivatives of this 
position), and its dynamics are specified by an equation 
of motion for (f ) . 

Consider now a quantum description of a superfiuid 
with a single vortex in it. One begins with the N- 
particle wave-function ^({r.,-}) — ({rj}^), where j = 



2 



1,2, ....TV and \rj) is a position state for the j-th parti- 
cle; or, equivalently, from the iV-particle density matrix 

PNiiTjhVi}) = ({Tj}\PN\Wi}), where p N = 
The state- vector |W) and the wave-function *&({]".,■}) are 
assumed to satisfy Bose symmetrization over permuta- 
tions of the particles. We now add a vortex soliton to 
the system, with circulation q v K, where q v = ±1,±2, 
and assume that in the iV-particle wave-function the 
vortex node is at a point R(t) in the plane. Both 
the iV-particle wave-function and the total density ma- 
trix pjv must then depend explicitly on the parame- 
ter R(t). One is now free to make a change of vari- 
ables to a set {R;q,t} of 'collective coordinates' [15] . 
with k = 1,2, ....N — 1, wherein the vortex coordinate 
R(t) is isolated from the remaining N — 1 coordinates 
{qfc}. This vortex coordinate, formerly just a parame- 
ter in the wave-function, is now elevated to the status 
of a quantum variable associated with a state vector |R) 
and with a 'zero mode', which we discuss below. All 
other degrees of freedom must now be properly orthog- 
onalized, both to the zero mode and to each other [T5] , 
The corresponding A^-particle density matrix is written 

Piv({R > qfc};{R , ,q' k }) = ({R;_q,-}|piv|{R , ,q , J -}). We 
also define a vortex reduced density matrix by integrating 
out the {q_k}- 

p v (R,R',t) = Tr qfc p w ({R,q fc };{R',q fc }) 

= II / dq k p N ({R,<lk};{R',<lk}) (1) 

ft 

in terms of which all physical quantities associated with 
the vortex may be defined and evaluated, provided we 
have sufficient information about Pn- 

One may now, by fairly well-established manoeuvres, 
derive a field theory for the Bose-condensed system, 
starting from the 1-particle reduced density matrix [13], 
denoted by pi(r, r', t). Note that this is a different object 
from the vortex reduced density matrix we have defined 
above, and is obtained in the usual way by integrating 
over the coordinates of all the particles in the full den- 
sity matrix except for one of them (in a fully symmetrized 
way that takes account of the Bose statistics). As is well 
known, one can characterize the result in terms of a quan- 
tum phase field $(r, t) and a density field p(r,t); they 
are assumed to have Bose commutation relations (and in 
a classical approximation, these fields are amalgamated 
into a macroscopic wave-function ip(r,t) ~ p(r, t)e l ^^ r,t \ 
and the commutation relations are dropped). It is com- 
mon in the literature to separate out a slowly-varying 
'texture' in the two fields, and write 



$(r,<) = $ s + 0(r,t) 
p(r,t) = p + rj(r,t) 



(2) 
(3) 



where the 'quasiparticle' variables (f>(r, t) and r](r, t) de- 
scribe fluctuations about the texture. However once we 



introduce a vortex into the superfluid, we have to be a 
little more careful. The vortex solution breaks the global 
translational symmetry of the superfluid, leading to a 
new quantum mode associated with this broken symme- 
try, the so-called vortex zero mode. Without specifying a 
particular Hamiltonian H for the superfluid, we can nev- 
ertheless say that the general equations that admit the 
vortex are: 
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for a fixed vortex. The perturbed Hamiltonian in the 
presence of the vortex is quadratic in the phase and den- 
sity variations <f>, rj and leads to the coupled equations: 
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A trivial solution to the perturbed Hamiltonian can be 
found by taking the gradient of the original vortex equa- 
tions Q: 
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(6) 
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where $v and pv are the 'texture' solutions in the pres- 
ence of the vortex. Comparing with the equations of mo- 
tion resulting from the perturbed Hamiltonian, we see 
that the derivative of the original vortex profile satis- 
fies them at zero frequency. The zero modes are then 
0o = V$v ■ n, r)o = V pv ■ n where the derivatives have 
been projected onto an arbitrary direction n. 

The zero mode can also be found by considering a small 
translation of the vortex. Expanding the vortex solution 
about a shifted position r + Sr, we have: 



$ v (r + Sr) w $y(r) + Sr ■ V<J>y(r) 



(8) 



= $ v (r) - Srsm(6 - 9 d )-d $ v (r) 
r 

pv (r + St) ps p v (r) + Sr ■ V p v (r) (9) 
= p v (r) + Srcos(9 - 8 d )d r p v (r) 

The zero mode corresponds to the prefactor of the small 
translation Sr, ie., the zero mode generates translations of 
the vortex, and indeed, the zero mode degrees of freedom 
are equivalent to the vortex degrees of freedom. 

There are two key observations that follow from this 
discussion. First, and rather obviously, we must now 
exclude the zero modes when defining the quasiparticle 
excitations. Thus, the quasiparticle wave-functions must 
now be redefined so as to be orthogonal at all times to the 
zero modes, so that as the vortex moves, the quasiparticle 
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wave-functions must continuously adapt to the changing 
position of the vortex (indeed, they are excluded from 
the vortex core, and phase-shifted far from the vortex). 
Second, although there will be an interaction between 
the vortex and the new perturbed quasiparticles, this in- 
teraction cannot have any term linear in the quasiparti- 
cle variables. This is because the vortex soliton is itself 
a minimum action solution to the equations of motion, 
and so any fluctuations about this solution (correspond- 
ing to the perturbed quasiparticles) are at lowest order 
quadratic in the fluctuation variables. 

These points are familiar in the discussion of quan- 
tum solitons in 1 + 1-d field theories [13]; a well-known 
example is the quantum Sine-Gordon model. However 
1 + 1-dimensional field theories are in many ways rather 
unique, and the standard belief for a long time has been 
that quantum soliton problems in higher dimensions were 
intractable. In fact this is not the caseJTOj; however, the 
vortex problem does bring in some interesting new fea- 
tures, notably: 

(i) unlike most of the interesting 1 + 1-dimensional 
models, the quasiparticle spectrum is gapless here. This, 
along with the long-range nature of the interaction be- 
tween vortices and quasiparticles, emphasizes the infra- 
red part of their coupling - indeed we expect to find diver- 
gences in the coupling to the unperturbed quasiparticles 
(which would mean that any perturbative or diagram- 
matic expansion in powers of this coupling would be at 
best unreliable, at worst meaningless). However, as we 
shall see, the coupling to the perturbed quasiparticles is 
not IR divergent. 

(ii) the perturbed quasiparticles differ from the unper- 
turbed plane wave quasiparticles not only in the spatial 
form of their wave-function - they are also now chiral ex- 
citations, with angular momentum defined relative to the 
vortex position. 

The upshot of all of this is that we must now distin- 
guish between the original quasiparticles, described by 
the field variables 4>(r,t) and r?(r,i), and the new per- 
turbed quasiparticles. We can describe the low-energy 
dynamics of the system by defining these new variables 
as excitations about the vortex texture, ie., we write: 

*(r, t) =$ v (t- R(i)) + 4>(r, t\R(t)) (10) 
p(r, t) = p v (r - R(t)) + fj(r, t\R(t)) (11) 

where the notation makes clear that the quasiparticles 
are tied to the vortex position. We can then write <j) and 
f) in cylindrical components centered at the instantaneous 
position R(i) of the vortex, which we write as 

4>(r, 9, t) = 4> lk (r) sm(uj k t + 19) (12) 
fj( r , 9, t) = -fjik(r) cos(w k t + 19) 

We see that the time and angular dependence cannot be 
separated in the perturbed quasiparticles; they are now 
chiral modes tied to the background vortex. 



From this discussion one might imagine that we can 
now completely forget about the original plane wave 
quasiparticles. In an isolated system this would indeed be 
the case. However in experiments one can do something 
rather interesting, which is to inject 'external' plane wave 
quasiparticles - in effect, one can irradiate the vortex with 
an external quasiparticle wind. These quasiparticles are 
not orthogonal to the vortex 'zero mode' wave-function, 
and they will interact linearly with it. Below, we discuss 
the experimental implications of this point. 



VORTEX-QUASIPARTICLE INTERACTION 

Formally, we may now write the expansion of the su- 
perfluid action in terms of the perturbed (tilded) quasi- 
particles in the form: 

S = S° v [R(t)} + SMM + AS^Mfj}] (13) 

where the zero mode is accounted for in R(<); here 
5°[R(t)] is the vortex action and S qp [{4>, fj}] is the 
quasiparticle action (both written in terms of perturbed 
quasiparticles), and ASl nt [{</>, f]}} is the interaction term 
(where the superscript indicates that it is quadratic in 
the perturbed quasiparticle variables). Rather than give 

(2) ~ 

a lengthy discussion of how AS 1 ^ [{(/>, fj}] is calculated, 
let us instead discuss the result, which can be portrayed 
in terms of the Feynman diagrams for the final form of 
the vortex-quasiparticle interaction. One may give these 
results either in terms of unperturbed 'external' quasipar- 
ticles, or in terms of the perturbed (tilded) quasiparticles 
- here we focus on the perturbed quasiparticles. 

The key question is of course to understand the form 
of the interaction between the vortex and the perturbed 

(2) ~ 

quasiparticles, which is incorporated in ASl nt {{<p, fj}. To 
give an intuitive feel for this interaction, we sketch here 
the form of the effective field theory which describes it, 
in diagrammatic terms. 

In a large system, the quasiparticle propagators are the 
same for perturbed or unperturbed quasiparticles. We 
define the quasiparticle matrix propagator G^ m (uj) by 

(hp s k 2 ninuj \ 
»w 4 > = 1 (14) 

In the same way one may define a propagator for the 
vortex itself, starting from 5°[R(t)]. Consider now the 
diagrams for the interaction between the quantum zero 
mode and the perturbed quasiparticles. One of the ver- 
tices involved is shown in Fig. 1. 

The expression for the total coupling Ajj£, between the 
vortex and a pair of quasiparticles having momenta k and 



4 



k v 



the quasiparticles in a Bose superfluid, and it is useful to 



P - (k+q) 




FIG. 1: Interaction between a vortex zero mode and a pair 
of perturbed (tilded) quasiparticles in a Bose superfluid. The 
dashed and wavy lines indicate quasiparticles of opposite chi- 
rality, and the momenta of each excitation are also shown. 



q respectively, takes the following form 
dr 



1V kq 



2m 



r(<j>lkd r fjl+<r,q + Vlkdr<t>l+<r,q) (15) 
+ a{l + cr) (</>/£ ?); +0%( j + f)lk<t>l+cr,q) 



This interaction is zero unless <j — ±1: the moving vor- 
tex thus only couples the renormalized modes 0;k and 
fji'q to each other if \l — l'\ = 1, transferring angular 
momentum ha = h{V — I) between them. Because of the 
long-range vortex field, we focus on the 'long-wavelength' 
regime where ka D < 1. In this regime, only the transi- 
tions between I — and I = ±1 contribute to linear order 
in a k. Using the anti-symmetry of A under exchange of 
the initial state (k,l) and the final state (q, Z + er), viz., 
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, we can fully express the total coupling 



by the I = term, to get: 



1V kq — 



5(k - q) 



(lY+aVk-q 



q(k-q) 



k+q 



a Q q 



k(k-q) 



(16) 

if k < q 
if q < k 



valid for q v = I. 

The form of this interaction is interesting. As noted 
above, it is not IR divergent - this is because the per- 
turbed quasiparticle wave-functions have adjusted to the 
vortex zero mode wave-function. Nevertheless it is not 
analytic about the zero momentum point, and part of it 
changes sign with the angular momentum transfer a. In- 
sofar as one believes the long-wavelength description of 
the Bose superfluid, the result is exact in the long wave- 
length limit. 

There have of course been many attempts in the past 
to derive the form of the interaction between a vortex and 



compare the result ( 16 ) with some of the forms found 
in previous work on this problem. These fall mainly 
into 2 categories. The first set of calculations attempts 
to calculate a scattering amplitude for quasiparticles in- 
teracting with a classical vortex potential. Such calcu- 
lations automatically yield a quadratic coupling, ie., a 
coupling to quasiparticle pairs. However, the long-range 
nature of the vortex field creates infra-red divergences 
in this scattering amplitude, which require careful dis- 
cussion [HI [TS]; moreover, the potential also carries an 
effective Aharonov-Bohm flux It actually turns 
out to be quite difficult to compare the details of such 
calculations with the results given here, mainly because 
(i) almost all of these calculations deal with the scatter- 
ing of plane wave excitations off a static vortex (with no 
recoil) ; and (ii) the vortex itself is not treated quantum- 
mcchanically. 

A second class of calculations employs a Hamiltonian 
of form UM 



H 



1 



2M, 



■[P-<7„A(r] s 



(17) 



( c kqk • R + ^ [Pk + ™kw£q£] 



where the vector potential A yields a 'field' V x 
A(r) = Trhp s z, which is responsible for the Magnus force. 
This Hamiltonian takes the Feynman- Vernon/ Caldeira- 
Leggett form [T7], with couplings Ck to quasiparticle co- 
ordinates {qk} (having conjugate momenta {pk}), which 
are linear in the {qk}- As noted above, such a linear 
interaction does exist between the vortex and the unper- 
turbed plane-wave quasiparticles. However, as we have 
already explained, no linear interaction to correctly or- 
thogonalized quasiparticles can exist. 

Note however that this does not stop us from writing 



down a Hamiltonian like (17) containing a linear inter 



action between the vortex and pairs of quasiparticles - 
such forms have been employed in other cases involving 
quasiparticle-soliton interactions. However the interac- 
tion now depends on the momenta of both quasiparticles, 
and is usually strongly temperature-dependent; consider- 
able care is needed to evaluate it. In our view the only 
reliable way to carry out such a manoeuvre is to first 
derive the interaction between the vortex and the true 
orthogonalized quasiparticles, as above, and then from 
this derive the form of the interaction to a bath of effec- 
tive oscillators. 

As shown in ref.|10). one can in fact carry though a 
fully non-perturbative derivation of the time dynamics 
of the vortex system, by dealing directly with the super- 
fluid action, which incorporates the infra-red convergent 



interaction ( 16 ) between the vortex and the perturbed 
quasiparticle field. By then integrating out the quasi- 
particles one finds an equation of motion for the vortex 
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reduced density matrix - there is no need to deal directly 
with the vortex scattering problem at all. 



EQUATIONS OF MOTION 

As we noted at the beginning of this paper, the classical 
description of a vortex involves an equation of motion for 
the classical coordinate H v (t) . As we will discuss below, 
the correct classical vortex equation of motion turns out 
to be 



MvTLv — fM — fqp — fac(*) — fff\t) 



(18) 



where f ac (i) is some driving force, M v is the vortex mass, 
fhi = PsK x (R-u — v s) is the Magnus force for a vortex 
with circulation k = zh/m, and the quasiparticle force 

fqp i s 



fqp = D Q (v n - R„) + D' Q z x (v n - R v ) 



(19) 



in which the longitudinal drag D a (T) and the transverse 
term D' Q (T) depend strongly on the temperature T. The 
classic discussion of Iordanskii yields 



D' Q {T) = —Kp n (T) 



(20) 



Finally, fif' (t) is a the classical limit of a 'fluctua- 
tional noise' term Ffi(t), whose behaviour is defined by 
its correlator ^ (i - 1', T) = (Fj^t, T)F j fl {t' , T)) . In the 
classical regime this correlator takes the form 



Xij(t-t',T) 



xg°(t- 



■f,T)~xl(r)M(*-0 ( 21 ) 



ie., it is entirely longitudinal, and displays Markovian 
white noise - we discuss the temperature dependence be- 
low. 



Now cqtn. (18) is in fact the original HVI equation^, 
[3] , but with an added inertial mass term and a longitudi- 
nal noise term. It is actually a local (in spacetime) equa- 
tion, ie., it can be written in the form C(t)Tl v (t) = f(t), 
where C(t) is a local differential operator acting at time 
£, involving forces and an inertial term which act on 
Ru(t) at time t only. There is actually no reason why 
the classical dynamics need to be local - one could eas- 
ily have, eg., 'memory' terms in the dissipation of form 
J dt'Tijit — t')R{(t') : and quite generally one could have 
an equation of form C(t, i')Ru(t'), where C is now some 
integrodifferential operator function of t and t' . How- 
ever, we will see that the correlation times in the classi- 
cal regime are very short, and a local equation is a very 
accurate approximation to the truth. 

Consider now the quantum dynamics. This has to be 
written in terms of an equation of motion for the reduced 
density matrix p„(R, R',t). Now the range of possible 
different forms for an equation of motion for Pu(R, R', t) 



is very large. Quite generally we might expect the equa- 
tion to be non-local in the variables R, R', and t. More- 
over, there is no reason to assume that we will be able 
to write the equation of motion in terms of simple forces, 
which are a classical notion of limited applicability in 
quantum mechanics. 

It is then refreshing to find that the actual time dy- 
namics of p u (R, R',t) do assume a fairly simple form, 
even in the quantum regime. For a discussion of experi- 
ments it is convenient to transform the equations of mo- 
tion to the frequency domain, and we discuss this in the 
next section. But it is also of interest to look at them 
in the real time domain, which we do here. The deriva- 
tion of these results is described in refs. [101 H2] • The key 
assumption is that we can make a Born-Oppenheimer ex- 
pansion, assuming the vortex velocity is small compared 
to the sound velocity in the superfluid. For the deriva- 
tion of the forces acting on the vortex this expansion is 
perfectly well behaved, and we can thus have confidence 
in the results. However the derivation of the effective 
mass is more subtle and the Born-Oppenheimer expan- 
sion misses the 'radiation reaction' terms (which also ex- 
ist classically). Consequently the Born-Oppenheimer ex- 
pansion yields a hydrodynamic mass M° for the vortex, 
without frequency-dependent corrections or higher time 
derivatives (eg., terms proportional to Rj). 

The results for the vortex dynamics can be written in 
terms of an equation of motion for p„(R, R',i), but it 
is more illuminating here to give them in terms of the 
dynamics of the arguments of (R, R'. We define the sum 
and difference variables R„ = g(R + R') and £ = R — 
R', so that Rt,(i) denotes a 'centre of mass' coordinate for 
the vortex, and £(t) a 'quantum fluctuation' coordinate 
about the centre of mass coordinate. 

One can then write the results in terms of equations 
of motion for these 2 variables. For the centre of mass 
coordinate one gets an equation which can be written as 



M°R„(i) - f M (Rv) 



Fq P [R„ 



F fluc (t) (22) 



where /j\/(R u ) is again the Magnus force, and where the 
new quasiparticle force Fq P is now a non-local functional 
of the velocity of the vortex, relative to the normal ve- 
locity. The equation for £(t) takes a somewhat similar 
form: 



M^(t)-f M (i)-F% P [i(t)}=0 



(23) 



where however now /a/(£) = PsQvK X £{t) (ie., this force 
is like the Magnus force acting on the relative velocity 
(Rt,(i) — v s ), except it now acts simply on the 'fluctua- 
tion velocity' £); where the 'quasiparticle' term does not 
depend on v„; and where there is no noise fluctuation 
term. 

We do not have space here to look in detail at the equa- 
tion of motion for but it is physically illuminating to 
look at the new quasiparticle force term Fgp[R„ — v n ]. 



() 



To make the connection with the HVI classical force /„. 



appearing in ( 19 ), let us write its quantum generalization 
as 



Ff[R„-v n ] 



v„] (24) 



where the idea is to separate out the two terms that in 
the classical limit lead to the drag force and the Iordanski 
force. 

Let us consider in detail the 'parallel' term F^Ru = 
v n ]. It is a functional of the prior vortex velocity parallel 
to the normal fluid; in fact it takes the microsopic form: 

Ff [R„ - v„] = A J2 (A^) 2 n kq (n k - n q ) (25) 

rnakq 

x I ds (R,(s) - v n ) cos[Q kq (t - s)] 

where we have assumed a quasi-2d film of thickness L Zl 
and the frequency fl k q is just the difference in energies 
between the two quasiparticles that interact with the vor- 
tex, ie., VL kq = Lu k — LtJq, where we expect that in the long 
wavelength regime, hui k = c s |k|, where c s is the sound 
velocity. We see that this force has just assumed a sim- 
ple 'memory' form, and that its instantaneous value and 
direction depend on the previous path traced out by the 
vortex (more precisely, the component of that the vortex 
velocity along that path that was parallel to v„). How- 
ever there is no requirement for F^ [R t , — v n ] at a given 
time t to be parallel to v„ at the position R„ of the vor- 
tex - if v„ varies with time or with position, then this 
will not in general be the case. 

Consider now the behaviour in time of this 'memory 
term'. To do this, let us imagine a vortex following a 
straight line trajectory - we can then write the parallel 
force as 



FiflR, 



d Sj D||(i-s;T)(v n -R,(s)) (26) 



where the function D» (t — s; T) is a non-local (in time) 
generalization of the HVI coefficient D a (T). Now this 
function turns out to depend only on the product kT(t — 
s) jh (so that its Fourier transform, as advertised, de- 
pends only on the ratio il — Ml/kT). In Fig. 2 we show 
its behaviour as a function of renormalized time. 

Suppose we write the longitudinal damping function 
D\\{t — s;T) as a sum of a local 5-function contribution 
and a retarded term. Now at first glance, one's intuitive 
expectation is that the 5-function term will simply be 
equal to the classical D Q (T) contribution. But this is not 
correct. Actually one finds that 

D {l (t-s,T) = ^-D (T)5(t-s) + 5D {l (t-s;T) (27) 

ii 26 11 

where the long-time tail term SD»(t — s;T) behaves 
roughly as ~ (t — s)~ 2 . Thus as one tends to low temper- 
atures or short times, such that fcsT(i — s)/h <C 1, one is 




1 1.5 2 

k B T(t - s) /n 

FIG. 2: The longitudinal damping time integral kernel Dm (t — 
s), shown as a function of the dimensionless variable ksT(t — 
s)/h: the coefficient decays roughly as (t — s)~ 2 at long times, 
a relatively slow decay coming from the linear low-frequency 
behaviour of D\\(Q). The arrow at the origin denotes a local 
in time (5-function) damping contribution. 



left with the 'quantum contribution' to the longitudinal 
damping - but this is 16 times smaller than D a (T), and 
moreover, it behaves like a local term! It is only at high 
temperatures or long times, when one integrates over all 
of -Dy (t — s,T), that the full contribution to the classical 
coefficient D (T) is recovered. 

From these remarks one sees that the crossover be- 
tween the classical HVI equations and the fully quantum 
regime is going to be an interesting one. We have no 
space here to outline the behaviour of the other terms 
in the equation of motion - suffice it to say that both 
the transverse quasiparticle force and the fluctuation 
force have non-local memory terms (although the cor- 
rection to the Iordanskii force in the transverse term 
turns out to be very small). Remarkably, once we have 
made the full crossover to the quantum regime (ie., where 
kT(t — s)/h -C 1, or where HQ <C kT, one actually ends 
up again with a local equation of motion for the vortex, 
this time coming only from the (^-function terms in the 
various memory kernels. This equation is 



M„R„ 



qp 



F„ c (t) 

(Q) 



r CQ) 

flue 



(t,T) (28) 



where the quasiparticle force Fqp is given by 



QP 



16 



R, 



■ D' a z x (v n - R„) (29) 



and where the fluctuation correlator is again Markovian 
and entirely longitudinal: 

X^\t- S ,T) = ^LxliT^Sit-s) (30) 

Thus the equation of motion in the extreme quantum 
regime has exactly the same form as the classical HVI 



7 



equation, but with quite different coefficients (except for 
the Magnus and Iordanskii terms, which have exactly the 
same coefficients). 

To summarize - one finds that the actual equations of 
motion for a quantum vortex are rather complicated ex- 
cept in the extreme classical regime, where they reduce 
to the standard HV1 equation, and in the extreme quan- 



tum regime, where they reduce to equation (28 1. The 



intervening crossover regime is expected to show quite 
different behaviour from either of the two limiting cases. 



CONCLUSIONS, AND REMARKS ON 
EXPERIMENTS 



experiments (where the bounce frequency is very high) 
(ii) experiments on the dynamics of single vortices in 2- 
dimensional 'pancake' cold BEC gases - this problem is 
discussed by Cox and Stamp [15]: and (iii) very low-T ex- 
periments on turbulence (where vortex motions can be 
extremely rapid). The detailed theory of such experi- 
ments remains an interesting challenge. 

This work was supported by funding from NSERC, 
from P1TP, and from CIFAR. The work benefited greatly 
from discussions at various times with (and encourage- 
ment from) David Thouless; we also thank Bill Unruh 
and Gordon Scmcnoff for useful comments. 



Condensed matter systems are populated by 3 differ- 
ent kinds of quantum excitation - extended quasiparticle 
modes, localized modes such as spins, or defects, and 
quantum solitons. ft is obviously of great importance to 
understand how these different excitations interact, and 
the debate over the nature of vortex-quasiparticle inter- 
actions has assumed a central importance in the theory of 
superfluids over the years. In this work and in refs.|10|,[T2] 
we present what we think is a solution to this problem, 
obtained by extending the theory beyond the purely clas- 
sical regime. We should however note some of the lim- 
itations of this work. First, it is only valid in the long 
wavelength limit - we ignore higher-order interquasiparti- 
cle interactions, and excitations like rotons in superfluid 
4He (in 4 He this confines us to T < 0.6 - 0.7 K). It 
also means that we cannot deal with the crossover in the 
vortex flow field between the 'near' regions where nor- 
mal fluid viscosity can be neglected, and the far region 
where the viscosity controls the flow. Thouless et al. [H] , 
have shown there are subtle problems involved in this 
crossover, and in fact we believe that the question of the 
total circulating normal fluid to be found around a vortex 
still needs to be settled. A second limitation is that the 
assumption of a slow vortex and a Born-Oppenheimer ex- 
pansion, so that we cannot capture all terms contributing 
to the inertial forces on a vortex. Finally, the work here 
describes vortices in a Bose superfluid - vortices in Fermi 
superfluids have to be dealt with separately. 

It is nevertheless interesting to speculate on how one 
might address the quantum regime in experiments. One 
needs high frequency vortex motion. The obvious candi- 
dates for experiments then include (i) vortex tunneling 
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